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1. Introduction 



The ALPHA collaboration has worked over the years on a determination of the QCD A- 
parameter starting from experimental low energy hadronic input and using perturbation theory 
in a renormalized coupling at sufficiently high energy scales. At these scales it was demon- 
strated that perturbation theory is very accurate. The quoted results for the MS A-parameter are 
A^ro = 0.60(5) [1] in the quenched approximation and A^ro = 0.62(4)(4) [2] with Nf = 2 
dynamical quarks. On the other hand the Nf = 5 value extracted by matching various exper- 
imental data to perturbation theory in the (not always very) high energy region translates into 
A^ro ~ 0.55 [3], when = 0.5 fm [4] is used. Superficially this suggests a nice agreement, but 

the perturbative matching across the quark thresholds [5] yields -^^/-^^ « 1.4 which does not 
connect smoothly to the Nf = 0,2 numbers. In order to say more about this comparison, the low 
energy scale ro should be replaced by an experimental observable and the continuum limit should 
be evaluated with a better confidence than it was possible in [2]. Significant progress in the under- 
standing of the continuum limit requires to simulate smaller lattice spacings with good accuracy. 
We will motivate this further in Sect. 2. The difficult simulations are the ones in large volume 
where for example the Kaon decay constant is to be determined to set the energy scale in GeV. We 
will briefly explain in Sect. 3 that our previous approach of using Schrodinger functional bound- 
ary conditions in that part of the calculation meets somewhat unexpected (practical) difficulties. 
Since these are, related to true dynamical fermion effects, they are theoretically interesting, but it 
appears to be better to switch to (anti)-periodic boundary conditions in this part of the calculation. 
In Sect. 4 we will finally discuss a determination of the dependence of the lattice spacing on the 
bare coupling. This represents a useful piece of information for fixing the parameters of the large 
volume simulations. 

Before entering our discussion we add a comment on the motivation. One might object to 
the whole project of a determination of the A-prameter that very precise lattice determinations 
for aj^(Mz) have already been published [6]. However, apart from the use of rooted staggered 
fermions, in these determinations perturbation theory has been used at rather low renormalization 
scales and for non-universal quantities (small Wilson loops). These are defined at the scale of 
the (lattice)-cutoff. It is apparent from the discussion in [6] that the use of perturbation theory is 
problematic. The known terms in the expansion either have large coefficients or, if one resums by 
choosing a different scheme, the renormalization scale becomes even smaller and the expansion 
parameter larger. In order to describe the data several higher order terms in the expansion are 
fitted. Thus it appears that a computation following the ALPHA-strategy, where the continuum 
limit is taken and perturbation theory is verified to apply for the considered renonnalized coupling, 
remains very well motivated. We do not see any alternative to this strategy if a full control of all 
systematics is desired. 

In the following considerations we use the standard 0(a)-improved theory with Wilson's gauge 
action and the non-perturbatively determined [7] coefficient Csw of the Sheikholeslami-Wohlert 
term [8]. 



2 



Nf = 2 simulations at small lattice spacings 



R. Sommer 



2. Cutoff effects in Za 

In [9] we have presented evidence that cutoff effects tend to be larger in full QCD than they 
are in the quenched approximation. Here we would like to draw the reader's attention to the non- 
perturbative determination of Za presented in [10]. It uses a Ward identity in the Schrodinger 
functional in an x 9/ 4L geometry with L 0.8 fm as in the quenched approximation [11]. It can 
be shown that the quark-propagator disconnected diagrams which enter the Ward idenity vanish 
in the continuum limit. They are of O(a^) at a finite lattice spacing. In contrast to the quenched 
approximation where already at a = 0.1 fm they were insignificant (in comparison to the numerical 
precision), for A^f = 2 they contribute an about 15% effect in Za at such a lattice spacing. Even if 
this effect disappears very quickly at smaller a, it is unpleasantly large at the lattice spacings one 
typically would like to include in a continuum extrapolation. 

In the mean time we have investigated the problem further, finding that qualitatively this effect 
persists if one changes the angle 6 in the spatial fermion boundary condition. Alternatively we 
considered the Ward identity between static-light states in such a way that disconnected diagrams 
are absent. Unfortunately, even when using HYP discretizations for static quarks [12] the statistical 
errors in Za become relatively large at the smaller lattice spacings. Still, we confirmed that Za 
defined in this way is rather close to the definition with light-light states but disconnected diagrams 
dropped. 

In general, cutoff effects are expected to be more prominent in correlation functions (and for 
time separations) where excited state contributions are very important. We therefore investigate 
at present whether the approximate ground state projection of [13] suppresses the disconnected 
contribution to Za thus accelerating the continuum limit. Whether this attempt is successfull or 
not, these difficulties suggest that one most likely needs smaller a with dynamical fermions than 
in the quenched approximation. We now turn to another strong effect of dynamical fermions - one 
that is expected to persist in the continuum limit. 

3. The large- volume Schrodinger functional 

Apart from the non-perturbative evaluation of renormalization constants, the Schrodinger 
functional also proved to be advantageous for the computation of hadron masses and matrix el- 
ements such as Fk in the quenched approximation [14]. A time extent of T = 3fm allowed to 
clearly isolate ground state contributions. We have then attempted to compute the pseudoscalar 
masses and decay constants for Nf = 2 with an L-' x T Schrodinger functional, keeping L > 2 fm 
and T « 2.5 fm. Indeed, at a quark mass around the physical strange quark mass (k" = 0.1355), the 
effective mass of the pseudoscalar correlation functions /a,/p (see e.g. [15] for their definition) 
exhibit short plateaux. An example is shown in the upper part of Fig. 1. 

However, the plateaux disappear quickly when the quark mass is lowered. For a quark mass 
of about half the strange quark mass {k = 0.13605), excited state contaminations are strongly 
present in both the vacuum channel and in the pion channel. The former yield contributions oc 
exp(— (r — xo)£'™'^) and the latter oc exp(— xq {Ef — m-j^)). Once these two leading contaminations 
ai^e included, fits to the correlation functions are still reasonable. We show a fit where we have 
fixed £'™'^ = Iniji, Ef = 3mji. These are the energies of the multi pion states with the coiTcct 
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Figure 1: The effective mass for the Schrodinger functional correlation function /p at j3 = 5.3 on a 24^ X 32 
lattice for K = 0.1355 and K ~ 0.13605. The fit described in the text is extended outside the fit-range as a 
dotted curve. The dashed line indicates the fitted pion mass. 



quantum numbers, when the interaction of the pions is neglected. At sufficiently large L this is a 
good approximation. 

We may conclude that multi-pion states are observed, as expected in the full theory. Their 
amplitude appears to be significantly stronger than with (point-to-point correlators and) periodic 
boundary conditions [16]. The standard Schrodinger functional boundary operators have a strong 
overlap with these states. Even though it is interesting to observe these strong effects of dynamical 
fermions and a consistent description over a significant range of xq can be achieved in the form of a 
fit, their presence hampers a reliable estimation of the systematic errors. We have hence decided to 
switch to periodic boundary conditions for the purpose of computing large volume matrix elements. 

4. The lattice spacing as a function of the bare coupling 

As a first step towards such computations we now compute, in a massless renormalization 
scheme, the dependence a'(^o) of the lattice spacing on the bare coupling go for 0.04 fm < a < 0. 1 fm. 
Of course the function a(go) is not unique, but only defined up to cutoff effects, which depend on 
the renormalized quantity that is held fixed. We employ a renormalization condition which is rela- 
tively easily evaluated and which does not introduce artificially large a-effects. This has proven to 
be the case for the standard Schrodinger functional coupling g^{L), defined in [17, 18], at vanishing 
quark mass. 

We further specify a scale L* by 

/(L*)=5.5, (4.1) 
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10 
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Table 1: Raw simulation results and interpolated values. Values of am 
estimated to be at most 5 x 10^^. 



indicate that |z| = \Lm\ is 



which is known to lead to L*/<3>8 for the planned range of a. For such a choice, table 7 of [2] 
shows a change of by about Af^ = 0.3 when the boundary 0(a) improvement coefficient Ct 
is changed from its 1-loop to its 2-loop approximation. Using the non-perturbative beta-function 
of [2], 

d 



L—g' = -Igm) = 0.21(l)f^ at g' « 5.5 , 



dL 



(4.2) 



a value Ag^ = 0.3 translates into a 5% change in L* and thus a. The definition eq. (4. 1) is completed 
by an exact definition of the massless point. We choose the PCAC mass m (with non-perturbative ca 
[13]) with Schrodinger functional boundary conditions, with T = L = L* , 6 = 0.5 and a vanishing 
background field. 

Good guesses for the bare parameters go, fc at a prescribed L/a are easily made starting from 
table 11 of [2]. When the result of a determination of g^{L) is close to the target eq. (4.1) and m 
is close to zero, we may correct by a first order Taylor expansion with derivatives eq. (4.2) and an 
estimate of 



1^ 
L dm 



g \l- 



From the results at two different values of m and fixed jS 



^/gl 



(4.3) 

5.3574 in Table 1 we extract 



i3 


log(LVfl) 


5.3000 


2.056(08) 


5.3574 


2.120(11) 


5.5000 


2.368(14) 


5.6215 


2.474(14) 


5.8097 


2.776(19) 


Table 2: 


Results for 


L*/a. 





5 = 2.2(5) at f{L)^ 5.5. 



(4.4) 



The rest of the simulation results of that table are then corrected to match 
the target with this value of s (including its error) and with eq. (4.2). We 
aiTive at Table 2 where a precision between 0.8% and 1.9% is seen. These 
numerical values are very well described by the simple linear interpola- 
tion formula 



log(L7a) =2.3338 + 1.4025 (i3-5.5) 
as seen in Fig. 2 where a ±0.02 "en^or band" is shown. 



(4.5) 
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Figure 2: The results for L* /a as a function of j3. 



5. Outlook 

Using the estimate a 0.08 fm at j8 =5.3 [16], we have estimated the pairs {^,L/a) = 
(5.5,32) and (5.7,48) in order to remain in the large volume region L > 1.9fm. We are currently 
carrying out first simulations at these parameters. Quark masses on the L/a = 4S lattice are initially 
designed to be only slightly below the mass of the strange quark. The reason is that our first goal is 
to carry out a precise scaling test, which is best done at not too small quark mass. Combining with 
the results of [16, 19] a significant range of a close to the continuum can be covered. 

The simulations are currently being done with the DD-HMC algorithm [20]. Release 1.0 of 
Martin Liischer's software [21] has been adapted for the BlueGene/L and an efficiency around 
30% has been achieved. The simulations do thus run at a sufficient speed to expect results from the 
BlueGene/L in Jiilich rather soon. These efforts are part of coordinated lattice simulations (CLS) 
carried out together with other lattice groups at CERN, Madrid, Mainz, Rome (Tor Vergata) and 
Valencia. 
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